Nuclear spin relaxation induced by a mechanical resonator 
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We report on measurements of the spin lifetime of nuclear spins strongly coupled to a microme- 
chanical cantilever as used in magnetic resonance force microscopy. We find that the rotating- frame 
correlation time of the statistical nuclear polarization is set by the magneto-mechanical noise origi- 
nating from the thermal motion of the cantilever. Evidence is based on the effect of three parameters: 
(1) the magnetic field gradient (the coupling strength), (2) the Rabi frequency of the spins (the tran- 
sition energy), and (3) the temperature of the low-frequency mechanical modes. Experimental results 
are compared to relaxation rates calculated from the spectral density of the magneto-mechanical 
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Sensitive detection of nuclear spin signals requires a 
sensor that can couple strongly to the weak nuclear mag- 
netic moment. For inductively detected nuclear magnetic 
resonance (NMR), where a resonant coil picks up the rf 
signal of the precessing spins, this is achieved by scaling 
down the coil size so as to improve the current-per-fiux 
ratio [1]. Optimized inductive detectors are currently 
able to observe ensembles containing roughly 10 12 proton 
spins, equivalent to about 10 8 net magnetic moments [2]. 
Stronger couplings enabling higher spin sensitivity have 
been shown with magnetic force sensors, such as used 
in magnetic resonance force microscopy (MRFM) [3, 4], 
which recently measured the net moment of about 10 3 
nuclei [5]. 

In order to detect spins with still higher sensitivity, 
the detector must be coupled even more tightly to the 
nuclear magnetic moment. For a mechanical detector 
sensitive to forces in the x-direction, this is achieved by 
increasing the magnetic field gradient which gener- 
ates the magnetic force F x = fi • where fj, is the 
total magnetic moment of the spin ensemble. Present 
MRFM technology employs gradients up to 10 6 T/m [5]. 
If the gradient can be successfully pushed into the range 
of 10 8 T/m, force detection of single nuclear spins may 
become feasible. 

The strong interaction between spins and sensor, how- 
ever, also increases the "back action" of the sensor on 
the spins and makes them more susceptible to detector 
noise. Here we are concerned with magnetic noise gen- 
erated by the thermally vibrating cantilever, but similar 
effects can be expected in any real- world detector, as for 
example in an inductively-coupled rf circuit. The interac- 
tion between a mechanical resonator and spins has been 
the subject of a number of theoretical studies, and is pre- 
dicted to lead to a host of intriguing effects. These range 
from shortening of spin lifetimes [6, 7], to spin align- 
ment by specific mechanical modes either at the Larmor 
frequency or in the rotating frame [8, 9], to resonant am- 
plification of mechanical oscillations [10]. 



In this letter, we report direct experimental evidence 
for accelerated nuclear spin relaxation induced by a sin- 
gle, low-frequency mechanical mode. The observed re- 
laxation has its origin in the random magnetic field 
B(t) = ^x(t) created by the thermal (Brownian) mo- 
tion x(t) of the cantilever in a large field gradient. The 
nuclear spin and mechanical oscillator degrees of freedom 
are well decoupled in strong static fields, since the spin 
precession frequency is orders of magnitude higher than 
typical cantilever frequencies. However, when an on- 
resonance transverse rf magnetic field B\ is present, kHz 
frequency noise that overlaps with the Rabi frequency of 
the spin will induce spin relaxation. The rate of nuclear 
spin transitions may be described by the rotating-frame 
relaxation time Ti p , 
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where u\ = is the Rabi frequency of the spin, 7 

is the gyromagnetic ratio, and Sb z {^i) is the (double- 
sided) power spectral density of the z-component of B(t) 
[6, 11]. 

For a transversely oscillating cantilever as shown in 
Fig. 1, Sbz{w\) is related to the spectrum of cantilever 

tip motion by Sbz{^i) = (ifjf) S x (ui), and will — if 
S x (cui) has thermal origin — be proportional to k^T. 
Hence, in a strong coupling regime where the spin lifetime 
is dominated by noise from the mechanical resonator, 
we expect an explicit dependence on the field gradient, 
the Rabi frequency, the temperature, and the mechanical 
mode spectrum. 

We find evidence for mechanically induced spin relax- 
ation while measuring nuclear spin correlation times for 
small ensembles of statistically polarized 19 F spins in a 
CaF2 single-crystal sample [12, 13]. The setup for these 
experiments, shown in Fig. 1(a), combines three compo- 
nents: (1) an ultra-sensitive cantilever, (2) a nanoscale 
ferromagnetic tip, and (3) a micron-scale rf circuit for 
rf field generation. Experiments are carried out with 
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FIG. 1: (a) Experimental arrangement (setup A). A CaF2 
single crystal containing 19 F nuclei is attached to the end of 
the cantilever and placed in a fixed position ^100 nm above 
a nanoscale FeCo magnetic tip. An rf magnetic field induces 
magnetic resonance in a thin "resonant slice" of spins where 
the Larmor resonance condition is fulfilled. Raising or low- 
ering the external field B ex t shifts the resonant slice position 
up or down, allowing spins further from or closer to the tip 
to be selectively addressed, respectively. 



two different arrangements: setup A uses a 90- /im-long 
single-crystal Si cantilever [14] together with a FeCo thin- 
film conical tip and an external rf microcoil [5]. Setup 
B employs a slightly longer cantilever (120 /im) com- 
bined with a FeCo cylindrical pillar integrated onto a 
lithographically patterned rf microwire [15]. The CaF2 
samples, a few /im 3 in size, are glued to the end of the 
cantilevers. 

For spin detection we rely on the "rf frequency sweep" 
method to drive adiabatic spin inversions, thereby modu- 
lating the z-component of the nuclear polarization at the 
fundamental cantilever frequency [15, 16]. Here, the rf 
center frequency u r f/27r is 114 MHz with peak frequency 
deviation of the rf sweep Au/2tt in the range of 400 to 
1000 kHz. We measure the correlation time r m of the 
cantilever tip oscillation amplitude, which reflects the 
correlation time of the nuclear polarization [13]. r m can 
be determined in several ways, for example by calculat- 
ing the autocorrelation function [13], by measuring the 
linewidth of the associated power spectral density, or by 
using a bank of filters of different noise bandwidths [4]. 
r m is very closely related to the rotating frame relaxation 
time Tip, as will be discussed below. 

In a first set of experiments we observe that spin re- 
laxation depends strongly on the distance between spins 
and magnetic tip. We attribute this observation to an 
associated variation of the magnetic field gradient. Fig- 
ure 2(a) plots the relaxation rate r" 1 as a function of 
applied external field £> ext . Since the external field sets 
the region of space where the resonance condition is met 




0.2 0.5 1 2 

magnetic field gradient G e \\ (G/nm) 



FIG. 2: (a) Spin correlation time r m as a function of external 
field Bext- Data points (dots) represent experimental values, 
measured with setup A. Spins are resonant at a total field of 
l^exti + Btipl = 2.85 T. (b) Effective field gradient G e s(B ex t) 
as a function of B eyLt calculated from tip model, (c) Relax- 
ation rate r^ 1 as a function of field gradient G e ff obtained by 
combining (a) and (b). Solid line is a best fit with the lowest 
point excluded. Dashed line is a guide to the eye. 



— indicated by the "resonant slice" in Fig. 1 — changing 
£>ext is equivalent to changing the region in the sample 
where we probe the spins. Spins located closer to the tip 
require less external field to satisfy the Larmor condition, 
^y\B ext z + £?ti P | = ^ r f- For these spins, the resonance ap- 
pears on the low end of the £> ext scan [Fig. 2(a)]. Like- 
wise, spins far away from the tip experience the smallest 
£?tip and require the highest £> ext . 

In order to relate r" 1 to the magnetic field gradient, 
we first calculate the tip field B t i P (r) as a function of 
position r. A good model for the tip can be obtained 
by inferring the tip size from a scanning electron micro- 
graph, and combining it with MRFM data for the tip 
moment [5]. To estimate the effective lateral field gra- 
dient G e ff(£>ext) oc for each resonant slice, we then 
derive the gradient from B t [ p (r) and average it over the 
resonant slice volume [17]. As a result, we can relate the 
gradient to the external field and use this knowledge to 
plot 7~ m as a function of G e ff, shown in Fig. 2(b) and 
(c). 

We find that spin relaxation increases with the gradient 
as r" 1 oc G^ 3±0 ' 16 . Since the gradient is the main pa- 
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FIG. 4: Power spectral density of the thermomechanical noise 
at T = 4.5 K, showing the first three vibrational modes of 
cantilever B. Sketches illustrate the bending of the cantilever 
beam. Frequency, spring constant and quality factor of the 
lowest mode are fi = 2.57kHz, ki = 86/iN/m and Q± = 
17, 000. 



FIG. 3: (a) Spin relaxation rate r" 1 as a function of rf field 
magnitude uj\ = 7B1, measured with setup B. Each time coi 
falls below a cantilever mode, indicated by dotted vertical 
lines, a new dissipation channel is added. At low fields r m is 
presumably limited by the adiabatic condition (dashed line, 
Ref. [21]). (b) Spin relaxation rate r" 1 as a function of mode 
temperature T mo de, here for the third cantilever mode. The 
thermal motion of the mode is found to affect spin relaxation 
when cji/27t < f$ (upper figure), but not when coi/2ty > fa 
(lower figure). Solid line is a best fit. 



rameter describing the coupling between spins and oscil- 
lator, a stronger gradient is equivalent to a stronger cou- 
pling — hence the same mechanical noise creates more 
magnetic noise. While gradient-induced spin relaxation 
is the dominant relaxation mechanism in Fig. 2, other 
(intrinsic) processes will eventually dominate in the limit 
of small gradient. In Fig. 2, this may be the case for 
gradients G e ff < 0.3G/nm. 

In a second set of experiments we investigate the spin 
relaxation rate as a function of rf field magnitude u\ = 
jBi. Here we are able to measure r" 1 for uji/2tt between 
60 and 170 kHz. The rf field magnitude is calibrated by 
a spin nutation experiment [15]. For small the spin 
correlation times are so short that the measurement is 
dominated by the response time of the actively damped 
cantilever (~ 20 ms). The upper limit on uo\ is set by 
heating restrictions of the rf microwire. 

The measured relaxation rates are shown in Fig. 3(a). 
We find that the spin relaxation rate decreases rapidly 
until reaching a distinct knee at uj\/2it w 125 kHz. We 
believe that this feature is directly related to the thermal 
vibration of the third cantilever mode, having a resonance 
frequency of fs = 122 kHz [see Fig. 4]. 

To further explore the role of the third mode we can 
modify its effective temperature, denoted by T mo ^ e . For 
that purpose we excite the mode with bandwidth-limited 



noise using a piezoelectric actuator mechanically coupled 
to the cantilever. At the same time as we measure r m , 
we also monitor the mean-squared motion x^ ms of the 
cantilever tip in that mode. We can assign an equiva- 
lent mode temperature T mo d e = &i£? ms /&B, where ki is 
the effective spring constant of mode i and /cb is Boltz- 
mann's constant. Note that while T mo d e can become 
very large for strong mechanical actuation, the actual (or 
bath) temperature of the cantilever remains at T = 4.5 K, 
as the heat capacity of the single mode is much smaller 
than the heat capacity of the cantilever's phonon bath 
and the modes are only very weakly coupled. 

We investigate the dependence of r" 1 on T mo d e at 
several different B\ fields and find the two different be- 
haviors shown in Fig. 3(b): If uj\/2tt lies above f%, no 
change in spin relaxation is observed even for strong ac- 
tuation. On the other hand, for uji/2it below ^3, relax- 
ation is greatly enhanced. We can fit the dependence of 
the relaxation rate on the mode temperature and find 
r" 1 oc (T mo d e )°' 43±0 ' 0T . Because the slope does not level 
off even for low T mo d e , the spin relaxation rate is set 
by the thermal fluctuations of the third cantilever mode 
over the entire investigated mode temperature range, in 
particular at thermal equilibrium where T mo d e = T. At 
uji/2it = 98 kHz, we can also try to enhance spin re- 
laxation by actuating the second mode. No influence is 
observed, as expected, because u)i/2tt > fi> 

In order to better understand the dependence of the 
spin correlation time on field gradient, rf field magnitude, 
and temperature, it is worthwhile to connect r m to the 
cantilever's mode spectrum. We calculate the transition 
rate following the analysis of Mozyrsky et al. [6]. To 
obtain an expression for r m similar to Eq. (1), taking 
into account the time-dependence of the rf field frequency 
uj T f(t) during cyclic spin inversion, we assume that r m 
can be described as the average relaxation rate over a 
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(See also Ref. [6], Eq. (7)). Here, T c = 1/ fa is the 
oscillation period of the fundamental cantilever mode, 
^eff(^) = {[cJrf(^) — ^o] 2 + uo\} 1 / 2 is the effective mag- 
netic field [18], and oo T t{t) is modulated from c^o — Acj 
to uoq + Acj. Note that because Acj ^> oj\, the Rabi fre- 
quency o; e ff traverses a broad range of frequencies during 
the sweep. A key difference between r m and T\ p is there- 
fore that r m is sensitive to noise in a frequency band set 
by u)\ < ou e ff(t) < Acj, while T\ p is influenced by noise in 

the vicinity of uj\ only. 

The magnetic noise spectrum Sbz(u) is given by 



S Bz M = G, 
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where n is the number of modes with a significant noise 
contribution. Because of the high quality factors Q^, 
Sbz(&) exhibits a discrete set of sharp peaks at the mode 
frequencies fa. We can facilitate the analysis of Eq. (2) 
by treating the peaks as J-functions, and find that 



7 2 G e 2 ff fc B T 
2Aw 



(4) 



where n' is the lowest mode whose resonance frequency 
f n t is above coi/2tt. In other words, when uo\ is reduced 
a new relaxation channel is added each time coi/2tt be- 
comes less than a cantilever resonance fa. While n can 
be large, the lowest mode n' will dominate (4) because 
it has the lowest effective spring constant and hence the 
largest thermal vibration amplitude. 

As a result, we find that Eq. (4) describes both the 
kink in Fig. 3(a) at the frequency of the third mode 
and produces the correct trends in G e ff and T mo de- Eq. 
(4) and Refs. [6, 19], however, predict that r" 1 oc 
Gg ff T mo de — a significantly stronger dependence on G e s 
and T mo de than our experimental observation, r" 1 ex 

G l ff 23±0.16 T a43±0.07 f We are aware Qf twQ possible rea _ 

sons for this discrepancy: First, the analysis of the spin 
transition rate [Eq. (2)] is based on the Bloch-Redfield 
approximation, where the magnetic noise is assumed to 
be uncorrelated [20]. This assumption may be violated 
because the correlation time of the noise, r c ^ = Qi/Trfa, 
can be significant, possibly even longer than r m . Further- 
more, we assume that spins are non-interacting, which 
will not necessarily be valid, especially at low rf fields 
where local fields can easily exceed B\. 

In conclusion, we find that nuclear spin relaxation can 
be induced by a single, low- frequency mode of a microme- 
chanical resonator. The spin relaxation rate is observed 
to increase both with the gradient and the effective mode 



temperature. We also find that long spin lifetimes are 
recovered when increasing the magnitude of the rf field 
to raise the Rabi frequency above the lowest mechanical 
mode frequencies. 
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